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Using a bound-state model of weakly bound gluons for glueballs made of two gluons and a natural 
generalization of the perturbative QCD formalism for exclusive hadronic processes, we present results 
for glueball production in radiative J/tp, T decays into several possible glueball states, including 
L / ones. We perform a detailed phenomenological analysis, presenting results for the more 
favored experimental candidates and for decay angular distributions. 
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I. INTRODUCTION 

The unambiguous experimental observation of two- or three-gluon bound states should be an important test of 
Quantum Chromodynamics and of its internal consistence. In the last years a huge experimental and theoretical 
work was in fact dedicated to this goal. Unfortunately, lower mass glueball states are supposed to fall in a mass 
range (a few GeV) which is already largely occupied by mesonic and/or hadronic states, including exotic states, 
hybrids etc. Moreover, glueball states are probably not narrow enough to allow a very clear observation. However, 
the intense experimental search [1, 2], complemented by lattice QCD calculations [3, 4] and several other theoretical 
contributions, like bag models [5], flux-tube models [6], QCD sum rules [7], has given evidence to a few possible 
candidates for two-gluon bound states. Typically, experimental searches for glueball candidates have been performed 
in favorable experimental environments. A good place to search for glueballs is in radiative J/ip, T decays, because 
of the "gluon rich environment" produced in the QQ decay, QQ — > 7<?<?, which should be particularly suitable for 
subsequent glueball formation. Two other processes of interest are central production in hadronic collisions and pp 
annihilation in the few GeV energy range. It is interesting to observe that one can take advantage also of negative- 
answer experiments, like e.g. hadron production in photon-photon collisions. In this case, meson production is 
expected to be favored with respect to glueball formation, which should only proceed via higher order (in the e.m. 
coupling constant) contributions. In some sense, the situation is reversed compared to radiative quarkonium decays. 
From a more quantitative point of view, this implies that the so-called "stickiness" parameter [8], S = [T(J/tl>, T — > 
7 72.) PS(jj — > 72)]/[r(77 — > 72) PS(J/tp,T — > 772)], where 72 is a hadronic resonance and PS denotes the phase 
space factor, should be much greater for glueballs as compared to ordinary mesonic states. The present status of 
the glueball search program can be summarized by recalling that there are two main candidates for glueball states 
[1, 2, 9]: a scalar, ++ state with mass M ++ ~ 1.5 or 1.7 GeV; a tensor, 2 ++ (or 4 ++ ) state with mass Mj++ ~ 2.22 
GeV; these observations are supported also by lattice QCD calculations [3, 4]. Other possible experimental candidates 
are [2]: a pseudoscalar, h state with mass M -+ = 2.14 GeV; a vector, 1 ++ state with mass M]++ = 2.34 GeV; 
a pseudotensor, 2 h state with mass M 2 -+ = 2.04 GeV; a 3 ++ state with mass M 3 ++ = 2.0 GeV. Present lattice 
calculations, however, predict larger masses for most of these states [3]. Other glueball candidates have also been 
proposed. In the following, we will not consider these additional states, limiting ourself to discuss some selected 
examples in detail; it will then be clear how to deal with the remainder. Additional information can be found in ref.s 
[1, 2] and references therein. 

Motivated by the fundamental role played by J ftp, T radiative decays in the field of glueball search, in this paper 
we study in detail these decay processes and analyze their implications for glueball production. 
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To this end, we use an approach previously developed by some of us and already applied to several processes 
involving the production and/or decay of hadronic states (including qq mesons and two-gluon bound states) at high 
energies [10-13]. This approach is a generalization of the formalism developed, starting from the early 80's, in the 
framework of perturbativc QCD for the study of exclusive hadronic processes at high transfer momentum [14]. Our 
model assumes that the dominant contribution to the production/decay of a hadronic state at large energy scale is 
given by its leading Fock state, considered as a weakly-bound state of valence constituents, with any internal angular 
orbital momentum. It has also been shown that this approach reduces to the usual perturbative QCD results in the 
case of L — states. Let us clarify from the start that, being based on perturbative QCD, this model is not suitable for 
predicting physical observables which are nonperturbative in their essence. Thus, for example, it cannot give sensible 
predictions for the masses and decay constants of light qq and two-gluon states. Rather, these quantities must be 
taken as input parameters of our model, obtained from experimental results or other nonperturbative approaches, 
like e.g. lattice QCD calculations or QCD sum rules. In some cases, like for decay-product angular distributions, it 
is possible to give predictions for physical observables which have little dependence on these nonperturbative inputs, 
as will be shown in the following. 

Let us remind that our approach was already applied to glueball production in radiative quarkonium decays in 
Ref.[10], in the case of a nonrelativistic distribution amplitude for the two-gluon bound state. Here we extend and 
generalize the analysis to the case of a generic glueball distribution amplitude; we also present analytical and numerical 
results obtained using several different models for the glueball distribution amplitude. 

We derive and discuss, in the framework of our model, a number of general analytical results for scattering am- 
plitudes, angular distributions etc., that may be useful independently of the numerical results presented here. This 
should hopefully facilitate the use of our approach for future, more detailed investigations along this lines. 

The paper is organized as follows: in section II we summarize our approach and present all relevant relations 
required. In section III we consider in detail the radiative decay J/ip, T — ► 7G, deriving the expressions of the helicity 
scattering amplitudes for several possible two-gluon bound states with S = 0, 1, 2, L = 0, 1, 2, 3, 4 and J = 0, 1, 2, 3, 4; 
we also derive the expression of the decay photon angular distribution in the quarkonium rest frame. Section IV is 
devoted to the presentation and discussion of several numerical results and predictions; we give our final comments and 
conclusions in section V. A few appendices contain most of the analytical results for the helicity scattering amplitudes 
and some useful analytical integrals required in the calculations. 

II. DESCRIPTION OF THE FORMALISM AND GENERAL RESULTS 
FOR HELICITY SCATTERING AMPLITUDES 

In this Section we outline the formalism utilized in the following calculations and present some general expressions 
for the scattering amplitudes. We limit ourselves to summarize the main steps of our approach. The interested reader 
will find more details in Ref.s [10-13], where the same approach was formulated and applied to several hadronic 
production/decay processes at high energy scales. 

The basic ingredients we need in order to evaluate observables for the process J/i/j, T — > 7G are the corresponding 
helicity scattering amplitudes, Ax^,x G -,x Q , where the A's indicate the helicities of the corresponding particles. From 
now on we always indicate by the subscript Q all quantities related to the generic QQ bound state, Q = J/ip, T. 

Notice that we know from first principles that the helicity amplitudes A\ y! \ G; \ Q can be written, in the quarkonium 
rest frame, as 

Ax^x a; x a (0 7 A) = K+a e tX ^ 4 e ,A 7 -A G CM . « 
where 7 , (f>j are respectively the polar and azimuthal angles specifying the photon outgoing direction in the quarko- 
nium rest frame, d 3 x x , are the usual Wigner rotation matrices and Ax 1 .x G are "reduced" amplitudes, depending only 
on A 7 , \q and on the dynamics of the decay process. Notice also that, since d 3 x x ,(8 — > 0) = <5a,a', Eq. (1) gives 

^A 7 ,A G ;A Q (0 7 = ^7 = °) = ^A T ,A G $\ Q ,\-,-\ G ■ (2) 

This means that we can evaluate the helicity amplitudes ^4a 7 ,A g ;A q in the simple kinematical configuration 9 y = 
(f) 1 = 0, find the reduced amplitudes Ax ,a g from Eq. (2) and utilize Eq. (1) to reconstruct the full amplitudes. 

In the framework of perturbative QCD and its factorization theorems for exclusive hadronic processes at high energy 
scales [14] (in our case, the QQ bound-state mass, Mq), the helicity amplitudes ^4a 7 ,a g; a q for the physical process 
may be evaluated starting from the helicity amplitudes for the elementary hard process QQ — > jg*g*, 7a 7i a 1 ,a 2 ;Aq,Aq , 
where Ai, A2 are the helicities of the two final (virtual) gluons. 
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The partonic process is schematically represented in Fig. 1, both in the QQ cm. ref. frame (the quarkonium rest 
frame) (A) and in the glueball rest frame (B), where the two constituent gluons are moving with relative momentum 
2k. As a first step, using the well-known nonrelativistic color-singlet model for the quarkonium bound state, we can 
easily find the connection between the helicity amplitudes Xa 7 .Ai.a 2 :A<3,Aq and those relative to the process Q — > 7 <?*<?* , 
■^A 7 ,Ai,A 2 ;Aq ■ A further step will in turn connect the amplitudes M to the hadronic amplitudes A, Eq. (1), which 
describe the full hadronic process. 

Let us remark at this point that in the last years the crucial role of color-octet contributions in hadronic processes 
involving the production/decay of charmonium bound states has been discussed in the framework of nonrelativistic 
QCD (NRQCD) [15]. It has been shown that color-octet contributions are crucial for the consistence of a nonrelativistic 
model of quarkonium. To our knowledge, only few attempts have been made to apply NRQCD and color-octet 
contributions to the treatment of exclusive charmonium decays [16]. These attempts have shown that color-octet 
contributions are probably of little relevance in the case of J ftp, T exclusive decays. In what follows, then, we stick 
to the usual color-singlet model for the 3 Si, QQ bound state. In this case we can write: 

M A7 , Al , A2;A J0 7 ^ 7 ;/3,M) = /% E C )£:t!?:L T ^MMM Q Wi,<l>*M,<t>), (3) 

where /q is the quarkonium decay constant; the C^' J M are the well-known Clebsch-Gordan coefficients; 6, (f> are 
the polar and azimuthal angles specifying the direction of the relative momentum between the two gluons, 2k, in 
the glueball rest frame; j3 is the modulus of the relative momentum in units of the glueball mass, j3 = 2\k\/Mc- 
Without giving more details on the calculation of the hard scattering amplitudes T, let us only recall that, at leading 
order in the strong coupling constant power expansion, there are six Feynman graphs contributing to the process, 
corresponding to all the possible ways the photon and the two gluons can attach to the fermionic line. The expressions 
of the helicity amplitudes Ma tj Ai,a 2 ;A q are collected in Appendix A. 

Following our approach [11, 12], the hadronic amplitude, A^ \ G J \ Q (#71 <P-y), Ior the production of a two-gluon bound 
state, G{gg), with spin, orbital, total angular momentum, and helicity S, L, J, Xq respectively, is given by 

^X, Q {e^) = 2 g l jjz^l}^ ECa^^^.a.a.^^,^;/?,^^), (4) 

Ai,a 2 

where c/l is the glueball decay constant, 

, _ 1 ( 2L+1 \ 1/2 (_ 2i \ L (2L + 1)!! [ d L 

Rl(t) is the glueball radial wave function in momentum space; <f>(z) is the glueball distribution amplitude, which 
describes how the glueball longitudinal momentum (in the so-called infinite momentum frame) is shared between its 
two constituent gluons (2:1,2 = (1± z)/2 is the fractional longitudinal momentum carried by gluon 1, 2 respectively), 
and 

& Xo (o, *) = <%fjc*x 8 i e " iAG * . (6) 

where A = Ai — A2. 

It can be easily checked that Eq. (4) reduces to the usual result of perturbative QCD for L = bound states, 
apart from some trivial redefinitions of the decay constants. In fact the same expression, with obvious appropriate 
modifications, can be used for a QQ bound state, and Eq. (3) could be retrieved from Eq. (4) by substituting the 
gluons with the QQ pair, putting L = 0, S = J = 1, and using for the J/ip distribution amplitude the nonrelativistic 
limit, & NR (z) = 6(z). 

Let us now consider in more detail the total glueball wave function, It must be completely symmetric under 
exchange of the two gluons, g\ < — ► 52- is factorized into the product of the color, spin, orbital wave functions, 
and of the distribution amplitude, = & c &s&l&(z)- Since the two gluons must be in a color singlet state, $ c 
is symmetric. The combination of two s = 1 particles can give S = 0,1,2, and $s(gi < — > 32) = (— l) s< i ) 5; the 
orbital angular momentum component has symmetry (— l) L ; finally, in our approach the distribution amplitude $(z) 
is always taken symmetric in z. As a consequence, to have a completely symmetric wave function, only states with 
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even L + S are allowed. We have explicitly checked that all amplitudes with odd L + S would require a nonvanishing, 
antisymmetric distribution amplitude. We do not consider this "exotic" possibility in the following. The parity of the 
glueball state is given by P = (— 1) L , while its charge conjugation is C — (— 1) L+S = +1 for all allowed states. 

A (pure two-gluon) glueball bound-state | J PC ) is given in the most general case by a linear combination of states 
with different values of L, S and same J, P, C quantum numbers (only states with even L + S and (— 1) L = P may 
contribute). In Table I we show, for each glueball state with quantum numbers J PC (up to J = 4) all the \L, S) states 
with L < 4 contributing. We can put this in a concise form: 



J+S 



fjrc\J PC )= E 7l + (-l) L+S l + P(~l) L A L , s f L \L,S), 



(7) 



S=0,l,2 L=\J-S\ 



where fjpc is the decay constant of the physical state and Al,s are (complex) mixing coefficients; whenever only one 
\L, S) state contributes, we take Al^s = 1 and fjpc = /z,. The expressions of the helicity amplitudes for all \L,S) 
states of interest are reported in Appendix B. 

Let us remark that in our a pproach the gluons in the two-gluon bound state have an effective constituent mass, 
mi,2 = (Mg/2)^/(1 ± z) 2 — /3 2 . This fact has two important consequences: i) The gluons have three possible helicity 
states, Ai,2 = ±1,0. It is easy to check (see Appendix B) that all amplitudes with Ai = and/or A2 = correctly 
vanish in the limit 7711,2 — * 0. ii) Yang's theorem states that two massless spin 1 particles cannot bind to form a J = 1 
state, or a state with odd J and negative parity. In principle this theorem does not apply to our model with constituent 
gluons; in fact we find nonvanishing amplitudes for the 1 ++ , 1 ^, 3 h bound states; however, all amplitudes involved 
correctly vanish in the limit m^2 — > 0. The observation of two-gluon, 1 ++ , 1 '", 3 h bound states in radiative decays 
of quarkonium should indeed be an interesting test in favor of approaches, like ours, involving massive constituent 
gluons [17]. 



III. DECAY RATES AND ANGULAR DISTRIBUTIONS 



In the previous section we have presented all general results regarding the helicity scattering amplitudes required 
for the calculation of decay rates, photon angular distributions, etc. In this section we first discuss the choice of 
the (nonperturbative) glueball distribution amplitude, $(z). We will consider two physically plausible and justified 
possibilities that, on one hand, give indications on the dependence of numerical results on <fr(z) and that, on the other 
hand, are such that the integrals in the expressions of the helicity amplitudes, see Appendix B, can be performed 
analytically. After this discussion, we will present the expressions of the decay rate and of the photon angular 
distribution for the process considered, in terms of the helicity amplitudes. 

Looking at the results of Appendix B, one sees that all the integrals appearing in the helicity amplitudes are of the 
following form: 



I' 1 $(z) z 2m 

where y — Mq/Mq. To perform numerical calculations, one needs an explicit expression for the glueball distribution 
amplitude, &(z). This function is largely unknown, given its nonperturbative origin. The simplest possible option 
would be that of considering the nonrelativistic limit, Q NR (z) = 5(z). This gives very simple analytical results for the 
helicity amplitudes; however, it is generally considered inaccurate, in particular for low mass bound states; therefore, 
we will not consider furthermore this option in the rest of the paper. Let us only notice that the reader can easily 
derive all useful results for the nonrelativistic case by himself from Appendix B. We will then consider two more 
realistic options, based on results obtained in the past years for the production/decay of light mesons at high energies 
and largely adopted in the literature in various forms: 

i) A generalized version of the so-called asymptotic distribution amplitude [13], 



*°">{z) =N L (1- zY +1 = ^( + X + 2) (1 " Z2)L+1 ' (9) 
where L(z) is the well-known Euler Gamma function. 

ii) Following ref.s[18-20], a version of the asymptotic distribution amplitude which suppresses end-point contributions, 
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<£ AS (z) = N u (l-z 2 ) cxp[- U /(l-z 2 )] 

= ^cxp(u/2) {u 2 K {u/2) + u{l-u)K 1 {u/2)Y 1 {I- z 2 ) cxv[-u/{l- z 2 )} , (10) 

where u = M G /(26 2 ), b is a hadronic scale parameter and K and K\ are the well-known modified Bessel functions 
of the second kind. Notice that when b — > (or better when Mg 3> b) & AS (z) — » <5(z); for light mesons, like e.g. the 
pion, and reasonable values of b <~ 0.5 — 1.0 GeV, $ j4s, (z) is very similar to the well known asymptotic distribution 
amplitude [14]. 

With these choices for the glueball distribution amplitude, all integrals ^i, m ,n{y), Eq. (8), can be performed ana- 
lytically. A more detailed discussion of these integrals and a collection of useful results are presented in Appendix 
C. 

The decay rate and the photon angular distribution for the process Q — > "fG{J PC ) can be easily obtained from the 
differential decay rate 



dT 



d</> 7 <2(cos0 7 ) 647r 2 M s 



(11) 



Aq , Xq A-y , Ag 



where p(Q) is the helicity density matrix of the quarkonium state, normalized so that Trp(Q) = 1. p(Q) depends on 
the quarkonium production dynamics. For typical production in high energy £ + £~ colliders all ofT-diagonal elements 
are vanishing, p±i,±i(Q) — 1/2, poo(Q) — 0; in the following applications we will consider this case. 
By using Eq.s (1), (2) one can easily see that the total decay width can be written in the form 



(12) 



while the decay polar angular distribution, normalized to unity, can be written as 



1 dT 



IMcos0 7 ) ^Ex y ,x G \Ax y ,x G \ 2 Xa ^ XG 



PA e ,Ae(Q) 



dia^-Xai -/) I-4a 7 ,a g | 2 - 



(13) 



If 1 4+i,o I ^ we can define, for a glueball state with total angular momentum J, 



J l^+i^l 2 



1 |i + i,o| 2 

where i = 1, 2 and, since |Ag| < J, rj = if J < i. Then, using Eq. (13), one finds 



(14) 



1 dT J 3 1- 2r{ + ri \ 1 



r J rf(cos6» 7 ) 4 1 



1 - 2r{ + r. 



7T3 PiAQ) 



i-3p M (e) 



cos 2 6L 



(15) 



If |4+i, o| = and |j4+i,i| ^Owe may in turn define 



\A 



+1,+2| 



(16) 



and Eq. (15) simplifies to 



dY" 



T J d(cos6» 7 ) 



4 1 + r-J, 



-Pi,i(e)- i-3pi,i(Q) 



cos 0~ 



(17) 
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IV. NUMERICAL RESULTS 



Since our approach cannot predict the glueball decay constants, /l, we are not able to give reliable estimates for 
absolute quantities like decay rates, unless using information from other nonperturbative approaches. Another source 
of uncertainty is that, as we have seen, even considering only pure two-gluon states, some of the bound states with 
given quantum numbers, J PC , may result from the mixing of several L, S states, with unknown mixing coefficients. 
Although for weakly bound systems it is usually assumed that amplitudes corresponding to larger values of L are 
suppressed, the overall contribution of each L, S state may be modified by the inclusion of the mixing coefficients. 
Therefore, in this paper we limit oursclf to give predictions separately for each L, S contribution and for quantities 
that do not depend on the glueball decay constants, i.e. the angular distribution of the decay products and the 
branching ratios R = T(T — > jG)/T(J/ijj — > jG). For this last quantity some experimental information is available 
concerning possible glueball candidates. Both these observables can in principle give useful information regarding the 
glueball distribution amplitude and, for a given J PC state, regarding the relative weight of the various L, S amplitudes 
contributing. 

A. Predictions for the ratio R = T(T -> 7 G)/T(J/tp — > 7 G) 

Let us consider the ratio R — T(T — > jG)/T(J/ip — ► 7G). i? is independent of the value of the glueball decay 
constant. On the other hand, it depends on the quarkonium decay constants fy, fa. From Eq. (12) and the results 
of Appendix B we easily see that 

„ = r(T^ 7 G) = l-y r ( a s (M*) \ 2 ff T \ 2 (M^\ 3 F(y r ) 

r(J/^ 7 G) 1-y, \a s (M*)) \u) \ My J F(y^) ' 1 ) 

where j/q = (Mg/Mq) 2 and F{y r ) includes the rest of the contribution from the helicity amplitudes not explicitly 
shown. We can make use of the well known leading order result 

rtfifft^e+O-p*^, (19) 
and of the leading order expression for the strong coupling constant 

127T 

" s(Q2)= (lln c -2n / )ln(QVA2) ' (20) 
with n c = 3, n/ = 4, 5 for Q = M^.My respectively, and A = 0.2 GeV. We then find 



R _ T(T^e+e-) ( \ 1 — y r ( 251n(M^/A 2 ) \ F{ y r) 

r(J/^eV) UtJ 1-2/* I 231n(M|/A 2 ) J F(yJ ' 1 > 



We make also use of the experimental results [9] 



r oxp (J/i/> -» e+e") = 5.26 ± 0.37 KeV, 

r exp (T^e+e-) = 1.32 ±0.04 ±0.03 KeV. (22) 

In Table II wc present our results for the ratio R, for several glueball bound states. For each state, the mass is 
chosen according to the most favored experimental candidate available [1, 2, 9]. Due to its interest, we also consider 
one of the states forbidden by Yang's theorem, i.e. the 1 ++ state, for which we have indicatively taken Mi++ = 2.34 
GcV. 

For each state, we have also separately considered all possible \L,S) states contributing. This could be useful, 
together with the study of the decay angular distributions, to get information on the unknown mixing coefficients. 
To study the dependence of our results on the choice of the glueball distribution amplitude, we present results for the 
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generalized asymptotic distribution amplitude ^ < £ AS , Eq. (9), and for the modified asymptotic distribution amplitude 
<I> AS , Eq. (10), using in this case two indicative values for the parameter b, b — 0.5 GeV and b = 1.0 GeV. 

For two of the most accredited glueball candidates [9], i.e. the resonances /o(1710) and /j(2220) ( J = 2 or 4), some 
experimental results on branching ratios are available: 



B cxp (J/^ -> 7 /o(1710)) x £? exp (/ (1710) - K+K~) = 8.5^ x HT 4 , 

B cxp (J/iP -» 7 /./(2220)) x B oxp (./j(2220) -» if+if ~) = (8.1 ± 3.0) x 10~ 5 , (23) 

B°*p(T 7 / (i7io)) x B oxp (/ (1710) -» < 2.6 x 10~ 4 , 

£ oxp (T -> 7 /j(2220)) x B oxp (./j(2220) -» #+#") < 1.5 x 10~ 5 . (24) 

Taking the lowest possible value for the J/tp case, this results in upper bounds for the ratio R, which are reported 
in Table II. 

From Table II we can see that, apart from two cases (i.e., the L = 2, S = contribution to the 2 ++ state and the 
L = 4, S = contribution to the 4 ++ state), the results show relatively little dependence on the glueball distribution 
amplitude, differing at most by a factor of two to three. Our estimates are consistent with the experimental upper 
limits, when available; for the 4 ++ state the results are very close to (and in some case slightly larger than) the 
upper limit, while in all other cases they are considerably smaller. For a given state, there are remarkable differences 
between the possible L, S contributions, which could be useful once more experimental information on the ratio R 
will hopefully be available. Complementary information can be found by looking at the photon angular distributions. 

Regarding the scalar, ++ glueball state, since the Jo (1500) is also reported to be a possible candidate as an 
alternative to the /o(1710), we have also considered the case M ++ = 1.5 GeV. The value of the ratio R correspondingly 
decreases (increases) by 10-20% (20-30%) at most for the L = (L = 2) contribution. 

In spite of the lack of detailed experimental information, we have investigated the dependence of the ratio R on 
the glueball mass for the other states as well, using other possible mass attributions [2]. While the h , 1 ++ , 2 ++ 
states show on the overall a dependence smaller than, or comparable to, that of the scalar case, the 2 h , 3 ++ , 4 ++ 
states show a stronger sensitivity. As an example, a 10% change in the mass of the 4 ++ state induces a change of 
up to a factor 2-3 in the ratio R, the sign and the value of the change depending on the particular L contribution 
and glueball wave function considered. This sensitivity could be of some usefulness in discriminating among possible 
candidates in the mass range where our model is applicable. 



B. Photon angular distributions 

Using the results of Eq.s (15), (17), we present in Fig.s 2-11 the normalized photon angular distribution, I(0 y ) = 
(1/r) dr/d(cos0 7 ), for the decays T, J/tp — > 7 G, considering the glueball candidates of Table II. For each J PC state 
we give the angular distribution for all possible L, S contributions separately, and for the three different choices of 
the glueball distribution amplitude. Notice that for the ++ and ^ states, where for each L, S contribution there is 
only one nonvanishing amplitude (the one with Xq = 0), I (6 7) = (3/8)(l + cos 2 # 7 ), which gives no useful information. 
There are also some cases in which, for a pure L, S contribution, I{6 1 ) is independent of the glueball distribution 
amplitude. This happens when, at fixed L, S values, all nonvanishing amplitudes with different values of Xg have the 
same functional form of the integrals over the variable z. More precisely, this situation applies to the following cases: 
the L — 0, S = 2 contribution to the 2 ++ state; the L = 1, S = 1 contribution to the 2 ^ state; the L — 2, S = 2 
contribution to the 3 ++ state; the L = 2, S = 2 contribution to the 4 ++ state; it also applies to the 4 h state, not 
considered in Table II. 

Regarding Fig.s 2-11 one can make the following remarks: 

1) Apart from the cases where, for the reasons discussed above, there is no dependence on the glueball distribution 
amplitude at all, practically all remaining cases show little or negligible dependence on &(z). The only remarkable 
case is the L = 3, S = 1 contribution to T — > 7G(2 '"), where there is a sizable difference between the results 
obtained with the generalized asymptotic DA and those with the modified asymptotic DA. 

2) Regarding the possibility of discriminating among the different \L, S) contributions to a given J PC state, the overall 
situation looks more interesting. There are remarkable differences in the 2 ++ case, for the J/tp decay; in the 2 
case, for the T decay; in the 3 ++ case, for both T and J/tp decays; in the 4 ++ case, for the J/tp decay. 

3) The 2 h , 3 ++ , 4 ++ cases are those that, on the overall, show the most interesting effects when comparing the T 
and J/tp decays. 
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V. CONCLUSIONS 

In this paper we have presented a detailed derivation of the helicity amplitudes, decay widths, and angular dis- 
tributions for two-gluon glueball production in radiative quarkonium decays. To this end, we have further extended 
an approach previously developed and applied by some of us to different high-energy hadronic processes. Given the 
limited knowledge of some of the nonperturbative ingredients entering the calculation (the glueball distribution ampli- 
tude and decay constant; the mixing coefficients among different L, S contributions) we limited oursclf to presenting 
some indicative numerical results for the decay widths and angular distributions in the production of possible glueball 
candidates in a pure L, S configuration, adopting two different models for the glueball distribution amplitude. These 
results are mainly intended to illustrate the potentiality of our approach in discriminating among different glueball 
candidates and their masses, quantum number attributions, etc., once it is complemented by experimental information 
and other theoretical, nonperturbative inputs. We have tried to present in detail all analytical results, so that they 
could be used and even further extended by the interested reader. 

An interesting feature of our model is that it implies the use of massive constituent gluons, so that two-gluon bound 
states can escape Yang's theorem and form 1 ++ , 1 '", 3 h states, whose observation in the mass range considered 
here would then be a strong argument in favor of approaches, like ours, involving constituent gluons. 

We conclude by stressing that our approach can be generalized to include light qq pair (mesonic) production and 
possibly exotic states made of a diquark-antidiquark pair. This could be a necessary step if one would consider more 
realistic situations, where the observed resonances may result from the mixing of qq, gg, etc. components. Let us also 
remark that a similar approach could be extended (and in fact this has been already done, at least in part) to the 
case of hadronic resonance production in photon-photon collisions. The completion of this unified approach would 
certainly be of great interest for a deeper understanding of hadronic structure and spectroscopy in the mass range of 
a few GeV/c 2 . 
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APPENDIX A: HELICITY AMPLITUDES FOR THE PROCESS J/ip,T -^19*9* 

In this Appendix we present the expressions of the helicity amplitudes, Ai,A 2 ;A c > f° r the process Q — ► 7 g*g* 7 
with the kinematics defined in the rest frame of Q, see Fig. 1(A). Notice that according to Eq.s (1), (2), it is sufficient 
to evaluate the helicity amplitudes in the simple kinematical configuration Q 1 — </> 7 = 0. Let us recall that, while 
the amplitudes are given in the quarkonium rest frame, the angles 9 and </> specify the direction of the two-gluon 
relative momentum in the glueball rest frame, see Fig. 1(B). A simple Lorentz boost along the z-axis connects the two 
reference frames. We further notice that, due to parity invariance, only 27 out of all 54 amplitudes are in principle 
independent. Therefore in the following relations we fix A 7 = +1. 

Let us define 

128 n^^ajQ 1 

9Mq [1 + y(f3 2 - z 2 )} 2 - (I - y) 2 ((3 cos 9 -z) 2 ' 1 ' 

where y — M^/Mq. We have for the helicity amplitudes Mx y! \ lt \ 2 -\ a (0, (j>) : 

M+, ±1±;+ = -k{{1 z 2 )(\ + cos 2 6) T 2y/3[(l T - z 2 } + 2z{2/3 ± y[(l T P) 2 ~ z 2 }} cosfl} , (A2) 

M +i±i=F; + = -K sin 2 8(1 +[3 2 - z 2 ) , (A3) 



M+,±,o;+ - T-V2 Ky/(1 - z) 2 - (3 2 sm9 (cos8± yz)(l + z) - (1 - y)f3 T V0 1 , (A4) 
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M + , ,± ;+ = ±V2K^(l + z) 2 - f3 2 sin0[(cos0 ± yz)(l - z) + (1 - y)/? ± y/3 



M +;0 ,o;+ = 2 *V(1 + ^) 2 - /?V(1 " " P 2 sin 2 , 



M + ,±, ±;0 = -V2X2/ 1 / 2 e-^sin^[(l - /3 2 - z 2 ) cos0 + z{2/?± y[(l =p /3) 2 - z 2 ]}} , 



M +i±iTi o = tV2 Ky 1/2 e _< * sin 0(1 =p cos6>)(l + /3 2 - z 2 ) 



(A5) 
(A6) 
(A7) 
(A8) 



M +i±i0;0 = T^y^V^-^-^e-^a T cos0) [(1 + z)(l + yz ± 2 cos0) T (1 - t/)/3 - y/3 2 l , (A9) 



M +;0 ,±;o - Mj^VM^e-^l ± cos9) \(1 - z)(l - yz =p 2 cos 6>) T (1 - y)/3 - y/3 2 l , (A10) 



M+,o,0;0 = -2x/2 J ftTy 1 /2 v / (1 + ^2 _ £2^ _ ^2 _ £2 e -** sinff cosff, 



(All) 



Af +t ±,± ; - = -Ay e _i2 * sin 2 0(1 - f3 2 - z 2 ) , 



(A12) 



M + ,±, T; _ = -Kye-WQ T C os0) 2 (l + /3 2 - z 2 ) 



(A13) 



Af+,±,o ; - = -V2Kyy/(l - zf - f3 2 e^iX + z) sin 5(1 T cos (9) , 



(A14) 



Af+,o,± ; - = -y/2Kyy/{i + zf - f3 2 e^il - z) sin (9(1 ± cos (9) , 



(A15) 



M +> o,o;- = -2lW(l + " /?V(1 " ^) 2 " /? 2 sin 2 . 



(A16) 



APPENDIX B: HELICITY AMPLITUDES FOR THE PROCESS Q — > 7 G(L, S, J) 

Using Eq.s (4), (5), (6), and the results of Appendix A, we here present the expressions of the helicity amplitudes 

^xl, S \g-\q (^7 — ^7 = 0)' f° r au ghieball states up to J = 4 and L < 4. The full angular dependence of the amplitudes 

can be obtained from Eq. (1). Moreover, due to parity invariance, A L J x s, f_ Xa ._ Xa = (_i)^-s+Aq-(A-,-a g ) A%' ^J. x . 

We here present only independent amplitudes, by fixing A 7 = +1. For each | J PC ) state we consider all possible \L,S) 
amplitudes contributing, according to Table I. 

JPC = ++ 



o.o.o 256V2 f Q \Rq(0)\ f 1 <Z> N (z) 1 

■A-l n 1 = T^T^V OtCt s — T77T I UZ — , ^ 

1,0,1 9V3 M Q M V2 vr^i-y 2 z 2 



(Bl) 
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JPC = -+ 



jPC = 1+ + 



2 . 2 .o 1024 /q |i%(0)| Z 1 , *n(z) 



9V3 M Q M V2 7_ 1 ^ T 372(l_ z 2)2 (1 _ y 2 z 2 ) 3 



x { 1 + I3y + y 2 + (3 + y - 8y 2 - 34y 3 + 3y 4 )z 2 

- y 2 (9 + 2y - 'Sly 2 - 5y 3 )z 4 + 5y 2 (2 - 3y)z 6 } . (B2) 



,1,1,0 -512 r- f Q \R'A0)\ f 1 , * N (z) 3-(l + 4y + y 2 )z 2 +3y 2 z 4 



2,2.1 2560V2 /g |i%(0)| ^ $n(z) 1 m ., 

2,2,1 _ 512V2 / a |fl 2 '(0)| 1/2 f 1 $n(z) 1 

1,li0 " 3V3 V S M S Mq' 2 ^ i-iV?(l-^(l-« 3 

x { 1 + 3y + y 2 - (2 - y - 7y 2 + 19y 3 - 3y 4 )z 2 + y 2 (6 - 17y + 16y 2 )z 4 } . 



. 512V2 /a |^(0)| i/ 2n , f 1 , 



. , VT^ 1 (l-z 2 )(l-y 2 z 2 ) 2 



jpc = 2 ++ 



,0,2,2 256 r- f Q \R (0)\ f 1 J 1 



(B5) 



.i, i, i . 5W2 /a |^(0)l M x f 1 , <Mf) 1 mfn 

A < * 1 = ? VT^Ma W (1 " y) li dz (i-z 2 )(i-^ 2 ) ' (B6) 



: 5i2^_ /a Ifl'iWLi^ „ A A 1 ^ i-y^ 2 (B7) 



A ' " 1 = -W^ as ^-M^ U z vi-. 2 1-^ ' (B8) 

4 o, 2, 2 256 / s |flo(0)| 1/2 f 1 <M*) 1 



0,2,2 256V2 /q |flo(Q)| f\ <Z> N (z) 1 

A, 2; -i = -— Vl _ ,2 ! _ ,2,. . 
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4>2>2 256V2 /g |flf v (0)| 1/2 



9V7 v S M Q M 9/2 y VT^(i-^ 2 ) 4 (i-y^ 2 ) 5 

x I -5 + 32 y + 108 y 2 - 40 y 3 - 32 y 4 

+ 2 (-l-86y-142y 2 + 518y 3 -1114y 4 + 796y 5 -160y 6 ) z 2 

+ (-25 - 84 y - 460 y 2 + 7232 y 3 - 16926 y 4 + 20156 y 5 - 11548 y 6 + 2760 y 7 - 160 y 8 ) z 4 
+ 4y 2 (31 + 311 y - 2408 y 2 + 4402 y 3 - 4199 y 4 + 1769 y 5 - 167 y 6 - 54 y 7 ) z 6 
+ y 4 (-686 + 3716 y - 4884 y 2 + 4256 y 3 - 1829 y 4 + 372 y 5 ) z 8 

+ 6y 6 (2- 130 y+ 107 y 2 -42y 3 ) z 10 + 63 y 8 z 12 | , (B18) 



,4,2.2 512 / a (Q)| r 1 $ w (z) 



'G 



9^7 V S M Q M 9/2 yr^2(l-z 2 ) 4 (l-y2 Z 2 )5 



J pc = 2" 



| 143 - 176 y + 120 y 2 - 32 y 3 + 8 y 4 

+ 2 (73 - 736 y + 1246 y 2 - 1256 y 3 + 700 y 4 - 256 y 5 + 40 y 6 ) z 2 

+ (31 - 144 y + 1792 y 2 - 2720 y 3 + 3618 y 4 - 2960 y 5 + 1768 y 6 - 480 y 7 + 40 y 8 ) z 4 

+ 4y 2 (-13 - 284y + 413y 2 - 784y 3 + 659y 4 - 380y 5 + 74y 6 ) z 6 

+ y 4 (362 + 208 y + 144 y 2 + 160 y 3 + 71 y 4 ) z 8 

+ 6 y 6 (-46 + 8 y - 25 y 2 ) z 10 + 63 y 8 z 12 } . (B19) 



lf !, 2 _ . 1024^ is_MM „ (1 _ dz Mg) (B20 ) 

Cx 

4 1, i, 2 . 512V2 /g |flU0)| 3 /2n ^2 [\ jMf) ^ 2 



x { 7 + (33 - 80y + 26y 2 )z 2 - (2 + 16y - 78y 2 + 16y 3 + 9y 4 )z 4 

+ y 2 (12 - 80y + 33y 2 )z 6 + 14y 4 z 8 } , (B23) 



, 3 ! 2 1024 _ fa \Ra'(0)\ i/ 2/1 ^ 2 f 1 , 



(B21) 



^;2;-i=0, (B22) 



.3,1,2 _ ,2048 r /q |i? 3 "(0)| n , 2 y 1 ^ <M*) i 



VT^2(l-z 2 )3(l-y 2 z 2 ) 4 
x 1 14(1 - y) + (28 - 157y + 188y 2 - 52y 3 )z 2 - y(3 - 144y + 142y 2 + 10y 3 - 18y 4 )z 4 

- y 3 (66 - 20y - lly 2 )z 6 + 21y 5 z 8 } , (B24) 



13 



3,i,2 ■ 28672 V2 /a |i%'(0)l n ,2 Z 1 , z 2 



jPC = 3 ++ 

A?;g;? = o, (B26) 

,2,2,3 1024 / a |i%(0)| 1/2 , 2 A 1 , *n(z) l+^^ + ^V+^V 

Cx 

2;2;3 102475 / a |i%(0)| f 1 <M*) l + (3-8y + 3y 2 )z 2 + y 2 z 4 

^ 2 - 1 = ^ 7r ^^^ y(1 " y) Lfy/l-z* (1-^(1-^ ' (B28) 

Cx 

,4, 2. 3 10240 r f Q \Rj v (0)\ 2 r 1 <M-*) l + (3-8y + 3y 2 )z 2 + y 2 z 4 
4,2.3 2560V2 / a |iff y (0)j 1/2 2 r 1 Z> N (z) 1 

= -^7T^ as M- Q -l^ v {l - y) i^^F^FiW 

x I 5 - 4y - 4y 2 + (20 - 40 y- 59 y 2 + 128 y 3 - 40 y 4 ) z 2 

+ (7 - 28 y - 152 y 2 + 664 y 3 - 737 y 4 + 260 y 5 - 20 y e ) z 4 
+ y 2 (-25 + 272 y - 716 y 2 + 632 y 3 - 169 y 4 ) z e 

+ y 4 (-43 + 140y-112y 2 + 24y 3 ) z 8 -3y 6 z 10 }, (B30) 

4,2,3 512V2 /a \Rj v (0)\ , 2 f 1 <t> N (z) 1 

K ^ = -—^ as M- Q ~l4r V{1 - y) Jj Z 7T=^(l-z 2 ni-y 2 z^ 

Cr 

x|53- 16y + 8y 2 + (200 - 712 y + 589 y 2 - 352 y 3 + 80 y 4 ) z 2 

+ (-5- 280y+ 1048y 2 -1136y 3 + 1063y 4 - 400y 5 + 40y 6 ) z 4 
+ y 2 (155 - 592 y + 664 y 2 - 712 y 3 + 215 y 4 ) z 6 

+ y 4 (-7+104y + 8y 2 ) z 8 - 15 y 6 z 10 }. (B31) 

JPC = 3 -+ 



3.4.3 . 512^7 /g I iff (0)1 2 r 1 



®n(z) 



VT^2 (l-Z 2 )3(l-y2 z 2)3 

{l - y + (1 - lly + lly 2 - 3y 3 )z 2 + (3 - y)y 2 z 4 } , (B32) 



3,1,3 . 4096^ /g |i%'(0)| 1/2 2 f 1 



®n(z) 



vr^(i-^ 2 ) 3 (i-r^ 2 ) 4 

{l - y - (1 + Uy - 16y 2 + 5y 3 )z 2 - y(3 - 24y + 20y 2 - 7y 3 )z 4 - y 3 (3 - y + 2y 2 )z 6 } , (B33) 
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3,1,3 . 8192V70 / a |i%'(0)| 2 f 1 $jy(£) 

^2-1=-' ™ 7iyaa s — 772-2/(1-3/) / <fe 7= i 



z 2 



i yr^(l-z 2 ) 3 (l-y^ 2 ) 



2 



(B34) 



jPC = 4 ++ 

2 ,2,4 1024V2 /affl,, , 2 *n(z) 1 + (3 - 8y + 3y 2 )z 2 + y 2 z 4 

^ i = ^2T^^i^ (1 - y) I^tct d-, 2 ) 2 (i-^ ■ (B35) 

2,2,4 1024V5 /a |i%(0)| 1/2 , 2 , gjKg) l + (3-8y + 3y 2 )z 2 + y 2 z 4 

i:i ' = ^r^M^i^^^ 1 -^ l^vi-z 2 d-z^d-^ ' (B36) 

2,2,4 1024V5 /a |i%(0)| . 2 Z" 1 . gg(f) 1 + (3 - 8y + 3y 2 )z 2 + y 2 z 4 

^^-^^T^Mal^^ 1 -^ (i-z 2 ) 2 (i-y 2 z 2 ) 3 ■ (B37) 
...... i 2V2 _ /q |i?f(o)| /- 1 ^ $ N (z) 1 



.4 



Aoi = F^Mt; —57^ / dz 



V3 ' M Q M 9/2 (l-Z 2 )4(l-y2 z 2 )5 

x { (1 - y) 2 (391 - 6102 y + 1721 y 2 ) 

+ (2503- 23466y+ 139087y 2 - 301024y 3 + 279356y 4 - 114996y 5 + 17210y 6 ) z 2 
+ y (-3442 + 67196 y - 413514 y 2 + 929047 y 3 - 915276 y 4 + 425324 y 5 - 89960 y 6 + 8605 y 7 ) z 4 
+ y 2 (625 - 40670 y + 320322 y 2 - 798058 y 3 + 773961 y 4 - 321024 y 5 + 52012 y 6 - 5788 y 7 ) z 6 
+ 2 y 4 (1955 - 46467 y + 153058 y 2 - 141519 y 3 + 42517 y 4 + 1096 y 5 ) z 8 

+ y 6 (20341 - 83242 y + 74871 y 2 - 23940 y 3 ) z 10 + 2660 y 8 z 12 | , (B38) 



^ 4 , , 4 _ 20480_ ^ /q |i? 4 /y (0)| 1/2/i A2 f 1 ^ 



^ /q 1^1" (0)1 1/2,, ,2 f , 
x{3-12y+12y 2 -4y 3 + (l - 20 y + 62 y 2 - 64y 3 + 28 y 4 - 4 y 5 ) z 2 



+ y 2 (6-20y+15y 2 -4y 3 ) z 4 + y 4 z 6 } , (B39) 



20480 /a |i? 4 Jy (0)| , 2 f 1 , jMz) z 2 [l + (3-8y + 3y 2 ) z 2 +y 2 z 4 ] 

4 ,2,4 5120 /a |fl 4 /y (Q)| 2 r 1 

^oji = --7==^^— /2 (1-y) / dz- 



(l-z 2 ) 4 (l-y 2 z 2 ) 3 
5120 ^ f Q |i?f(0)| ^ 2 y 1 1 



(B40) 



3VTT55 V M Q M 9/2 v »> J_, (l-z 2 ) 4 (l-y 2 z 2 )5 
x | 25 + 16 y - 8 y 2 + (96 - 176 y - 247 y 2 + 264 y 3 - 80 y 4 ) z 2 

+ (15 - 128 y - 144 y 2 + 1240 y 3 - 1101 y 4 + 400 y 5 - 40 y 6 ) z 4 
+ y 2 (-25 + 528 y - 1600 y 2 + 1408 y 3 - 597 y 4 + 88 y 5 ) z 6 

+ y 4 (-107 + 368y- 272y 2 + 88y 3 ) z 8 - 11 y 6 z 10 }, (B41) 



15 



4 , 2 , 4 512V2 f a \R{ V (Q)\ 1/2 2 f 1 $ N (z) l 

x I 3 (51 + 8y-4y 2 ) + (540 - 1452 y + 9 y 2 + 308 y 3 - 120 y 4 ) z 2 

+ (-5 - 500 y + 488 y 2 + 2564 y 3 - 1877 y 4 + 600 y 5 - 60 y 6 ) z 4 
+ y 2 (155 + 528 y - 3236 y 2 + 2508 y 3 - 1165 y 4 + 220 y 5 ) z 6 

+ y 4 (-287+ 1124y- 672y 2 + 220y 3 ) z 8 -55y 6 z 10 }, (B42) 
4,2,4 512V2 ,fg |fif v (0)| n , 2 f 1 , *n{z) 1 

A i 2-i = - - ^t W""^ — . .q/9 v\ l -v) I dz 



3V23T V S M Q M 9/2 ^ ^ 7-1 (1 - Z2)4(l _ y2 z 2)5 

x| 3 (79 - 48y + 24y 2 ) + (1160 - 4488 y + 5061 y 2 - 3168 y 3 + 720 y 4 ) z 2 

+ (195 - 3160 y + 10712 y 2 - 14064 y 3 + 10767 y 4 - 3600 y 5 + 360 y 6 ) z 4 
+ y 2 (995 - 4048 y + 5336 y 2 - 4488 y 3 + 1215 y 4 ) z 6 

+ y 4 (17 + 296 y + 72 y 2 ) z 8 - 55 y 6 z 10 | . (B43) 



jPC = 4 -+ 



3;1;4 .16384 /g l^( )l n , 4 /- 1 J <M*) z 2 [l + (l-4y + y 2 ) z 2 + y 2 z 4 ] 
^ = 8 ~ ^Ma^f^-^ L dZ 7T^ (1- Z2) 3 (1 _,2 Z 2 )4 . (B44) 



3,1,4 _ .20480V2 r f a K'(0)| 3/ 4 y 1 z 2 [1 + (1 - 4y + y 2 ) z 2 +y 2 z 4 ] 



3,1,4 .20480V2 / g fl;" 3/2 4 T $at . 



., (l-z2)3(l-y2 z 2)4 



(B45) 



•l : ;:J: ' = 0. (B46) 

APPENDIX C: USEFUL ANALYTICAL INTEGRALS 

In this Appendix wc collect several relations useful for the complete analytical calculation of the amplitudes 
^A 7 ,A G :A a - An inspection of the results presented in Appendix B shows that all the required integrals can be cast in 
the form of Eq. (8), 



*,„,„(„) - » jf * ^ „_X^ ■ (C1) 

We consider separately the two classes of glueball distribution amplitudes utilized in this paper. 



1. Generalized asymptotic distribution amplitude 

Using the generalized asymptotic distribution amplitude &1 AS (z), defined in Eq. (9), we easily see that the integrals 
*;,m,n(y) reduce to 

rl (1 _ 7 2\l/2 2m 

*f£%(y) = 2N L / dz [ — = N L l m ,a{y) , (C2) 

jo u — y z ) 



where I m ,n{y), independent of L, is given by 
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\ {l-z 2 )V 2 z 2m 



f 

z m , n {y) = 2 J dz ^ _ 



y 2 z 2y 

V^r(m + i/2) ,„ „ 2n 

= 2v\ m + 2) ;V i(m + l/Um + 2 lB »), 

and 2-^1(0, c, z) is the well known hypergeometrical function. It is easy to verify that 



(C3) 



Zmfi(y) =Zm,n(y = 0) 



7T (2m -1)!! 



2 m+l ( m + l)! 

Two useful recursive relations between the I m ,n{y) are the following: 



(C4) 



l m +l,n+l{y) = 2^^ 2 ' m ' n ^ ' n > °: 



(C5) 



1m+l,n+l{y) — -o Im,n+l(!/) ~ ^m,n(y) 

y l 

From these two relations one can easily get as well the following one: 



(C6) 



i m , n +i{y) = ( i + -Jr -7-) i m ,n{y) , n>0. 

2n ay / 



(C7) 



Eq. (C3) is not very useful for practical calculations. One can instead utilize the well known transformation formula 
for the hypergeometrical functions, 



T(c) T(-c') 

2 F 1 (a,b 1 c;x) = ±> [ J 2 F 1 (a, b, c 1 + 1; 1 - x) 
T{a')T{b') 

where a' = c — a, b' = c — b, and d = a + b — c, together with the so-called Gauss recursion formulas, 



(C8) 



2 F 1 (a,b,c;x) = 2 F 1 (a + l,b - l,c;x) + x 2 F x {a + 1, b, c+ 1; a;) , (C9) 

c 

2 F 1 (a,b,c;x) = ~ ^ ( ° ~ ^ ~ | j 2 Fx (a - 1, b, c - 1; a) + _ ^ _ 2 F i (a - 1, b + 1, c; x) , (CIO) 

in order to express I mjl (y) as a combination of several hypergeometrical functions of the simple form 2 Fi(r, d, d; x) = 
(1 — x)~ r , with arbitrary real d and integer r. 

For completeness, we present below all the functions T m ,n(y) required in the evaluation of the helicity scattering 
amplitudes of interest. 



2o,i(y) = 


7f , 

y 2 


Vi- 


y 2 ), 






(Cll) 






Vi- 


y 2 - 






(C12) 


= 


y 


Vi- 


y 2 _ 


1 2 

2 y - 


1 ^ 


(C13) 


isAy) = 




Vi- 


y 2 - 


1 2 

2 y - 

1 2 

"2 y - 


1 4 1 6\ 

8 y ~ 16 y ' ' 
8 y 16 y 128 y ; 


(C14) 


ZiAy) = 


-Vi- 


_ y 2 _ 


(C15) 



t n - - 4 ~ 3y2 

~ 8 (1-y 2 ) 3 / 2 
^i,3(y) = 



(1 -2/2)3/2 ' 

x 2 , 3 (y) = ^(l- 8(1 _ y2)3 /2 ). 



J 3,3(y) - 2^(6-1/ - 4(1-^)3/2 )■ 

t < \ 77 (ah 10 2 4 48 - 84j/ 2 + 35y\ 
^,3(2/) = ^To( 4 8-12y -2/ - (1 _ y 2 )3 /2 ) 

7T 8 - 12y 2 + 5y 4 

^0,4 (y) = 



^i,4(y) = 



16 (1 - 2/2)5/2 
7T 2-j/ 2 



16 (1 - t/ 2 ) 5 /2 ' 
X 2.4(y) = ^,2)5/2 > 

7T 16 - 402/ 2 + 302/ 4 - 52/ 6 

J 3,4(y) - ^( 16(1 _ 3,2)5/2 !)' 

7T 64 - I682/ 2 + 140i/ 4 - 352/ 6 2 
^.4(y) = ^( g(1 _ y2f/2 8 + y ). 



7T 64 - 1442/ 2 + 120i/ 4 - 352/ 6 



128 (1 - 2/ 2 ) 7 / 2 

7T 16 - 16?/ 2 + 5j/ 4 

128 (I-2/ 2 ) 7 / 2 ' 

7T 8 - 32/ 2 

128 (1 - 2/ 2 ) 7 / 2 ' 

7T 5 

128 (I-2/ 2 ) 7 / 2 ' 

7T 128 - 448?/ 2 + 5602/ 4 - 280j/ 6 + 35?/ 8 

£10 V 1 128(1 - ? / 2 ) 7 / 2 } ' 

* 2 640 - 2304j/ 2 + 30242/ 4 - I68O2/ 6 + 315y 8 . 
(10-» ~^tB )■ 



2y 12 v y 64(1 - y2)7/2 

7T (120 _ 20y 2_ y 4 1920 - 70402/ 2 + 95042/ 4 - 55442/ 6 + 11552/ 8 



82/ 14 v y y 16(1 - 2/ 2 ) 7 / 2 
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2. The modified asymptotic distribution amplitude 

Using the modified asymptotic distribution amplitude $ AS (z), defined in Eq. (10), the integrals ^L,m,n(y) take the 
form 

, ^ n ,T f 1 , exp[-u/(l- z 2 )]z 2m ,_ o , 
*«, m ,„(|/,«) = 2JV„ j o fe (1 _^,!l /2(i ;; ¥) „ . (C38) 

with the normalization factor 



3 

Let us now consider the integrals 



N u = - cxp(u/2) [u 2 K (u/2) + u(l - u)K 1 (u/2)] \ (C39) 



f 1 , exp[-u/(l - z 2 )] z 2m ,„ ln . 
J> > ( " ' = 2 i dz a _ z 2y% {1 _ y 2 z 2 )n ■ ( C4 °) 





It is easy to verify that the following recurrence relations hold: 

Jl, m +l,n+l(y,u) = ( Jl, m ,n+l(y,u) - Jl, m ,n{y, u) ) , 

Jl+l,m+l,n{y,u) = Jl+l,m,n{y,u) ~ Jl,m,n{y,u) , 

Jl,m+l,n+l(y,U) = — ^ -^Jl,m,n(y,u) , n>0, (C41) 

Jl+l,m,n(y,u) = —rJl,m.n{y, u) , 

aw 

Jl,m,n+l(y,u) = ^1 + ^ Jl,m,n(y^ U ) ' n>0. 

By defining i + 1 = 1/(1 — z 2 ) and 7 = 1/(1 — y 2 ), the integrals Ji, m ,n{y, u) can be cast in the form 



/■ 00 

Ji >ro ,„(j/,M)=e- u 7 " / d« m - 1 / 2 (t + l) i - m +"- 2 (t + 7 )-" e - tlt . (C42) 
Jo 

Solving these integrals is not immediate. However, one can show that, by performing a certain number of appropriate 
integrations and differentiations with respect to u (the order and the number of these operations being related to the 
values of I, m, n) and using as boundary conditions the simpler integrals J^i, m ,n(y,u = 0), they can be related to an 
integral representation of the well known error function, Erf (a;). 

Without entering into more details of this conceptually simple but algebraically cumbersome procedure, we present 
below a minimal set of three integrals of the family J^i^ m ,n(y,u) that, using the recurrence relations (C41), generate 
all other integrals required by the calculation of our scattering amplitudes. Since the expressions found are in some 
cases quite lengthy and involved, we do not present all of them explicitly. 

Jofi,o{y,u) = V^Ze-"-7r(u-l/2) [l-Erf(>/u)] , (C43) 
Jo,o,i(y,u) = J { 1 - Erf (Vu) -^l-y 2 cxp [ uy 2 /{\ - y 2 ) } [ 1 - Erf (y/u/{l - y 2 )) ] } , (C44) 

Jo,2,i(y,u) = ^e- tl {-2V^y 2 (4 + y 2 + 2V)+ [S+(8u-4)y 2 + (4u 2 +4u-l)y 4 ] ^ e 11 [ 1 - Erf (y/u) } 

- 8vMl-y 2 ) exp[«/(l - y 2 )] [ 1 - Erf [y/u/{l - y 2 )) ] } . (C45) 
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TABLE II: Numerical values of the ratio i? = T(T — > jG)/r(J/tp — ► 7G) for the glueball candidates considered in this paper; 
for each case, different \L, S) contributions (up to L = 4) and three possible choices for the glueball DA (see text) are considered; 
comparison with experimental bounds, when available, is also shown. 



\J PC ) 


M G (GeV) 


\LS) 




R 




^exp 








GAS 


ASI 














b = 0.5 GeV 


b = 1.0 GeV 




0++ 


1.71 





2.1 10~ 2 


2.2 10~ 2 


2.2 10~ 2 


< 3.4 10" 1 






2 2 


3.1 10~ 3 


2.1 10~ 3 


2.8 10~ 3 




cr+ 


2.14 


1 1 


3.8 10~ 4 


3.4 10~ 4 


3.5 10~ 4 




l++ 


2.34 


2 2 


2.8 10~ 4 


4.3 10~ 4 


3.6 10~ 4 




2++ 


2.22 


2 


7.7 10~ 3 


8.8 10~ 3 


8.5 10~ 3 


< 2.9 10" 1 






2 


1.3 10~ 2 


2.2 10 _1 


4.5 10~ 2 








2 2 


7.6 10~ 3 


6.7 10~ 3 


7.1 10~ 3 








4 2 


2.1 10~ 2 


7.7 10~ 3 


1.8 10~ 2 




2-+ 


2.04 


1 1 


1.3 10~ 3 


1.8 10 -3 


1.6 10" 3 








3 1 


1.5 10~ 2 


1.6 10~ 2 


1.7 10~ 2 




3++ 


2.00 


2 2 


1.4 10~ 2 


9.7 10~ 3 


1.2 10~ 2 








4 2 


1.0 10~ 2 


5.6 10~ 3 


1.0 10~ 2 




4++ 


2.22 


2 2 


3.2 10 _1 


2.010 -1 


2.710 -1 


< 2.9 10" 1 






4 


5.510 -1 


2.8 10~ 3 


1.410 -1 








4 2 


3.1 lO" 1 


1.5 10" 1 


2.8 10" 1 






FIG. 1: Kinematical configuration for the process QQ — > 7G in the quarkonium rest frame (A), and for the elementary part 
process QQ — > 7 g* g* in the glueball rest frame (B). 
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FIG. 2: Normalized photon angular distribution, 1(6-,) = (1/r) dT/d(cos # 7 ), for the T — » 7G(1 ++ ) decay and different 
choices of the glueball distribution amplitude (DA), asymptotic DA (ASY DA), b = 0.5 GeV (solid); ASY DA, b = 1.0 

GeV(dashed); generalized asymptotic (GAS) DA (dotted). 
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FIG. 3: Normalized photon angular distribution, 1(8-,) = (1/r) dT/d(cos # 7 ), for the J/tp — ► jG(l ++ ) decay and different 
choices of the glueball distribution amplitude (DA), $(z): asymptotic DA (ASY DA), b = 0.5 GeV (solid); ASY DA, b = 1.0 
GeV(dashed); generalized asymptotic (GAS) DA (dotted). 
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FIG. 4: Normalized photon angular distribution, 7(# 7 ) = (1/T) dr/d(cos# 7 ), for the T — > yG(2 ++ ) decay and different \L,S) 
glueball states, using the asymptotic DA (ASY DA), with b = 0.5 GeV: L = 0, S = 2 (solid); L = 2, S = (dashed); L = 2, 
S = 2 (dotted). Notice that the L = 4, S = 2 case (which is not shown, for clarity) is almost indistinguishable from the L = 2, 
S = case. Moreover, notice that for the L = 0, S = 2 case 7(0 7 ) is independent of the glueball DA; for all other cases, use of 
the ASY DA with b = 1.0 GeV, or of the GAS DA, leads to very similar results. 
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FIG. 5: Normalized photon angular distribution, 1(0-,) = (1/r) dT/d(cos # 7 ), for the J/tp — ► ^G{2 ++ ) decay and different 
\L, S) glueball states, using the asymptotic DA (ASY DA), with b = 0.5 GeV: L = 0, S = 2 (solid); L = 2, S = (dashed); 
L — 2, S — 2 (dotted); L = 4, S = 2 (dot-dashed). Notice that for the L — 0, S = 2 case I(# 7 ) is independent of the glueball 
DA; for all other cases, use of the ASY DA with 6 = 1.0 GeV, or of the GAS DA, leads to very similar results. 
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FIG. 6: Normalized photon angular distribution, = (1/r) dT/d(cos 7 ), for the T — > jG(2~ + ) decay and different \L,S) 

glueball states: L = 1, S = 1 (solid); L = 3, S = 1, ASY DA with b = 0.5 GeV (dashed), ASY DA with b = 1.0 GeV (dotted), 
GAS DA (dot-dashed). Notice that for the L = 1, S = 1 case 1(0^) is independent of the DA adopted. 
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FIG. 7: Normalized photon angular distribution, I{9 1 ) — (1/r) dT/d(cos 7 ), for the J/ip — > 7G(2 + ) decay and different 
\L,S) glueball states: L = 1, S = 1 (solid); L = 3, S = 1, ASY DA with b = 0.5 GeV (dashed), ASY DA with b = 1.0 GeV 
(dotted), GAS DA (dot-dashed). Notice that for the L = 1, S = 1 case I(6j) is independent of the DA adopted. 
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FIG. 8: Normalized photon angular distribution, 1(9-,) = (1/T) dF/d(cos 7 ), for the T — > yG(3 ++ ) decay and different \L,S) 
glueball states: L = 2, S = 2 (solid); L = 4, S = 2, ASY DA with b = 0.5 GeV (dashed), ASY DA with b = f .0 GeV (dotted), 
GAS DA (dot-dashed). Notice that for the L = 2, S = 2 case 1(0^) is independent of the DA adopted. 
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FIG. 9: Normalized photon angular distribution, 1(8-,) = (1/r) dr/d(cos# 7 ), for the J/tp — » 7G(3 ++ ) decay and different 
\L,S) glueball states, using the ASY DA with b = 0.5 GeV: L = 2, S = 2 (solid); L = 4, S = 2 (dashed). Notice that for the 
L = 2, S = 2 case 1(6-,) is independent of the DA adopted, while for the L — 4, S = 2 case the three choices here adopted give 
almost indistinguishable results. 
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FIG. 10: Normalized photon angular distribution, 1(9-,) = (1/T) d Y/d(cos 9-,), for the T — ► 7 G(4 ++ ) decay and different \L, S ) 
glueball states, using the asymptotic DA (ASY DA), with b = 0.5 GeV: L = 2, S = 2 (solid); L = 4, S = (dashed); L = 4, 
S = 2 (dotted). Use of the ASY DA with b = 1.0 GeV, or of the GAS DA, leads to very similar results. 
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FIG. 11: Normalized photon angular distribution, 1(9-,) = (1/T) dY / d(cos9-,), for the J/tp — ► jG(4 ++ ) decay and different 
\L,S) glueball states, using the asymptotic DA (ASY DA), with b = 0.5 GeV: L = 2, S = 2 (solid); L = 4, S = (dashed); 
L = 4, S = 2 (dotted). Use of the ASY DA with b = 1.0 GeV, or of the GAS DA, leads to very similar results. 



